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We introduce a density functional formalism to study the ground-state properties of strongly- 
correlated dipolar and ionic ultracold bosonic and fermionic gases, based on the self-consistent 
combination of the weak and the strong coupling limits. Contrary to conventional density functional 
approaches, our formalism does not require a previous calculation of the interacting homogeneous 
gas, and it is thus very suitable to treat systems with tunable long-range interactions. Due to its 
asymptotic exactness in the regime of strong correlation, the formalism works for systems in which 
standard mean-field theories fail. 


Introduction - In contrast with its widespread use and 
success in areas as diverse as quantum chemistry pQ , ma¬ 
terials science J2j or semiconductor nanostructures [3], 
Density Functional Theory (DFT) has received relatively 
little attention in the very active field of ultracold atomic 
gases. It is well known that the Hohenberg-Kohn theo¬ 
rems, originally formulated in terms of the electron gas 
mis], hold for both fermionic and bosonic systems, as 
well as for interactions different than the Coulomb one. 
However, the lack of adequate density functionals has 
hindered the role of DFT in the study of ultracold atomic 
gases in favour of other well-established approaches, such 
as the widely used Gross-Pitaevskii (GP) method in the 
case of Bose gases. The latter is a mean-field approach 
and does not allow treating the effect of correlations, 
which play a crucial role in many different phenomena 
occurring in ultracold quantum gases [BJ. One then often 
turns to configuration-interaction (Cl), quantum Monte 
Carlo (QMC) or Green’s-function methods (for recent re¬ 
views, see, e.g., Refs. EHSD- 

The advantages of DFT are very well known from 
electronic-structure calculations m- Being an in prin¬ 
ciple exact theory (although in practice relying on ap¬ 
proximations), DFT allows to go beyond the mean-field 
description by taking into account correlations between 
the interacting particles. The Kohn-Sham (KS) map¬ 
ping of the many-body problem into a non-interacting 
one makes it possible to apply the method to particle 
numbers orders or magnitude larger than those acces¬ 
sible with wave-function methods B SHE]: as well as 
to capture a large part of the effects of quantum statis¬ 
tics on the kinetic energy. Initial efforts have already 
been made to generalize the formalism to bosonic and 
fermionic ultracold quantum gases HM3- However, the 
biggest challenge in Kohn-Sham DFT -for both elec¬ 
tronic and ultracold atomic systems- is the construc¬ 
tion of good approximations for the so-called exchange- 
correlation functional [Tj, the term in the total energy 


describing the many-body effects beyond the Hartree 
level. The simplest approximations are those based on 
homogeneous interacting models, where analytical ex¬ 
pressions for the exchange-correlation energy per parti¬ 
cle are often available by fitting QMC mm or Bethe- 
Ansatz calculations na nn ng. Such so-called local- 
(spin-)density approximations (L(S)DA) have been gen¬ 
eralized and applied to the study of Bose mmm and 
Fermi naira ill ultracold gases with short-range inter¬ 
actions in different geometries and dimensionalities with 
promising results. The downside of L(S)DA-based ap¬ 
proximations, however, is that they become unreliable 
for systems in which the interactions between the parti¬ 
cles largely dominate over the kinetic energy and char¬ 
acteristic strong-correlation phenomena arise. Moreover, 
L(S)DA approaches may become unpractical for systems 
with tunable interactions, since for each different inter¬ 
action a previous many-body calculation of the uniform 
system with that same interaction is needed, which can 
be demanding even in the simplest one-dinrensional case 
[TS1 . The tunability of the interactions is relevant, for ex¬ 
ample, in the interesting case of ultracold dipolar quan¬ 
tum gases . These systems have recently received 

considerable attention, particularly for low dimensionali¬ 
ties, where the collisional instability towards head-to-tail 
alignment of the dipoles can be suppressed (see the re¬ 
views I22HM] ). Remarkably, dipolar interactions also al¬ 
low the control of reaction rates in ultracold chemistry 
(see Refs. [23 26]). One-dinrensional systems have been 
of interest in particular in connection with the discussion 
of crystalline structures in the strong coupling regime 
(see, for example, Refs. 

In this Letter, we introduce an alternative approximate 
functional to study ultracold gases with long-ranged in¬ 
teractions, based on the exact formulation of the strong 
coupling limit of the Hohenberg-Kohn density functional 
[33l 34j . This provides an effective single-particle po¬ 
tential in a rigorous and physical way [55M571 . without 
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relying on calculations of the uniform system energy by 
means of other many-body approaches. The formalism 
becomes asymptotically exact in the limits of both van¬ 
ishing and extremely strong coupling. The latter is obvi¬ 
ously the most interesting, because of the plethora of 
important phenomena that are out of reach of mean- 
field theories and problematic for wave-function meth¬ 
ods. Furthermore, our construction can be equally ap¬ 
plied to fermionic and bosonic gases, by simply changing 
the kinetic-energy functional. 

Formalism - The key idea of Kohn-Sham Density Func¬ 
tional Theory (KS DFT) is that the ground-state par¬ 
ticle density n(r) and energy E 0 of a general TV— body 
system with interparticle interaction Wj n t(r — r') and ex¬ 
ternal confining potential w ex t( r ) can be mapped (with 
some mathematical subtle caveats, see, e.g. (313) into 
a non-interacting problem with the same particle den¬ 
sity n(r), moving in an effective potential VKs(M> r ) = 
Vext (r) + I'Hxc ( H ; r) . The Hartree-exchange-correlation 
(Hxc) potential r'Hxc ([*i]; r ), which is a Lagrange mul¬ 
tiplier for the density constraint, is obtained from the 
functional derivative of the difference in energy 

Euxc [ n ] between the interacting and the non-interacting 
systems. Physically, Uhxc(M; r ) transforms the many- 
body interaction effects on the density into a single¬ 
particle potential. The non-interacting system is usually 
chosen in order to capture the relevant effects of statis¬ 
tics on the kinetic energy: for both fermions and bosons 
it is defined as the system with density n(r) and min¬ 
imum possible kinetic energy with fermionic or bosonic 
statistics, respectively. The problem is then reduced to 
the self-consistent solution of the Kohn-Sham equations 
[—|V 2 + v ext (r) + 'C H xc([™];r)] ( Mr) = Sifa(r), where the 
Hxc potential u Hxc (W;r) = i; mf ([n];r) + v xc ( [n] ; r) is 
the sum of the Hartree mean-field (mf) and exchange- 
correlation (xc) contributions, the latter needing to be 
approximated. The KS single-particle orbitals fa deter¬ 
mine the ground-state density of the system via the rela¬ 
tion [5] n(r) = JTrij \fa(r)\ 2 , where rq is the occupancy 
of the itli orbital. For bosonic systems at zero temper¬ 
ature one has no = TV, and by neglecting the exchange- 
correlation term in the KS potential one recovers the 
Gross-Pitaevskii equation, widely used for the study of 
dilute ultracold Bose gases with short-range interaction, 
where the effects of many-body correlations do not play 
an important role. 

Here we introduce the “strictly-correlated particles” 
(SCP) functional fd®{p p [n], which is complementary to 
the KS non-interacting kinetic energy: it is defined as 
the minimum possible expectation value of the particle- 
particle interaction in a given density n(r). For the 
Coulomb interaction, this functional has been widely 
studied [331 US HD], and it has been shown to be able 
to capture the physics of the strongly-correlated regime 
in model quantum wires and quantum dots [35M37j , yield¬ 


ing results beyond the mean-field level. The construction 
of for a given density n{ r) is equivalent to an 

optimal transport (or mass transportation theory, a well- 
established field of mathematics and economics) problem 
with cost given by the interaction mi an. While several 
rigorous results have appeared recently in the mathemat¬ 
ics literature I43H50| . here we provide a simplified physi¬ 
cal overview. The idea is that if we minimize the expec¬ 
tation of the interparticle interaction in a given density 
n(r), we must have a non-zero probability to find one 
particle wherever n(r) ^ 0. The many-particle state is 
then a continuum superposition of “strictly-correlated” 
configurations (iq = r,r 2 = f 2 (r),..., rjv = fjv(r)), with 
r spanning the whole region where n( r) 0, 

|'Fscp(ri,...,r Ar )| 2 = ^ / dr ^ d(ri - f p(1) (r)) 

P 

x <5(r 2 - f p( 2 )(r)) • • -S(r N - f p(JV) (r)), (1) 

and p denote permutations of {1,..., TV}. The co-motion 
functions {f^} are highly non-local functionals of the den¬ 
sity satisfying the equations mmm 

n(r)dr = n(fj(r))dfj(r) , (2) 

which are simply equivalent to state that the probability 
of finding particle “1” in the volume element dr around 
r is the same as that of finding particle “i” in the volume 
element dfj(r) around f,(r). The {f;} also satisfy cyclic 
group properties that ensure the indistinguishability of 
the particles mm- The SCP functional is then HDU 5T1 

1 r N 

^int CP N = 2 / drn(r)J2 v int(r - f*(r)). (3) 

Notice that here we explicitly consider the possibility 
of anisotropic interactions depending on r r' (with 
Uint(r) = Dnt(—r)) and not just on |r — r'|, such as the 
interaction between dipoles aligned by an external field 
MU. The same formal steps of Refs. [331 3S ID] can be 
repeated for this more general kind of interaction, leading 
to the following exact equation for the functional deriva¬ 
tive vscp([»];r) = ^£ (r) lKl 

N 

V« SC p(H;r) = V%t(r - f,:(r)), (4) 

i —2 

which has a clear physical meaning: the potential ^scp ( r ) 
represents a force field equal to the net interaction felt by 
a particle at position r due to the other TV — 1 particles. 
Our “KS-SCP DFT” approach consists in using [n] 
to approximate the mean-field plus exchange-correlation 
terms of the total energy functional or, equivalently, its 
functional derivative i>scp(M; r ) of Eq. @ to approxi¬ 
mate the Hartree-exchange-correlation potential 133133, 

Uks(M;t) ~ Uext(r) +uscp(N;r) . 


(5) 
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This way, both the kinetic energy and the many-body 
interactions are treated on the same footing in the self- 
consistent KS equations. 

A few remarks are necessary on the kind of interac¬ 
tions t’i n t(r) for which the KS-SCP DFT can be ap¬ 
plied. Several rigorous results are available for convex 
repulsive long-ranged interactions depending on |r| only 
[331 13 ESHSD]. general, for the SCP formalism to 
be physically useful, the interaction u; n t(r) needs to be 
long-ranged, otherwise the SCP solution of Eq. ([!]) is 
just one of the many minimizers (and actually the one 
with maximum kinetic energy) for the interaction alone 
in a given density (see Ref. [52] for a discussion on con¬ 
tact interactions). The SCP functional is thus naturally 
very well suited for ionic gases in the strong-correlation 
regime, where it is expected to provide a large part of 
the total interaction energy 13511371 . Even more interest¬ 
ing is the case of general dipolar anisotropic interactions: 
the SCP functional combined with the KS kinetic energy 
(fermionic or bosonic) should be able to capture many 
of the interesting phenomena observed in the strong- 
correlation regime [23 [33] . In this case, the SCP solution 
can be constructed from the dual Kantorovich formula¬ 
tion EH, for which few results have started to appear 
recently [531 55] ■ 

Applications to low-dimensional dipolar ultracold gases- 
We consider N ultracold bosonic or fermionic particles 
with dipole moment d in quasi-one- (Q1D) and quasi- 
two-dimensional (Q2D) geometries. We model these sys¬ 
tems with the external harmonic potential u ex t(r) = 
\(ijj 2 x 2 + u 2 y 2 -\-uj 2 z 2 ) in the cases where, respectively, 
uiy, oj z ui x , and w z oj x = uj v = w_l (effective 
Hartree units are used throughout the paper). For these 
geometries, assuming that all the dipoles are oriented 
in the same direction due to the action of some exter¬ 
nal field, one can derive effective Q1D and Q2D dipole- 
dipole interaction potentials (see e.g. Refs. [29] and 
[551] for the explicit expressions) by integrating out the 
harmonic motion along the very strongly confined di¬ 
rections from the three-dimensional potential Vdd { r ) = 
d 2 (l — 3 cos 2 d r d)/r 3 , where 0 T & is the angle between the 
dipole moment and the relative position between two 
particles. Since the anisotropic interaction requires the 
development of a dedicated dual Kantorovich algorithm 
that will be the object of future work, here for the Q2D 
dipolar systems we restrict ourselves to the case in which 
the dipoles are perpendicular to their plane of motion and 
their interaction is purely repulsive and isotropic [53 1. 

In order to perform practical calculations with the KS- 
SCP approach, one proceeds as follows. First, for a given 
density distribution, the co-motion functions { f,} are ob¬ 
tained by solving Eqs. ([2]) and performing an angular 
minimization [361 [37]. The one-body potential uscp( r ) 
can then be obtained by integrating Eq. <©• Finally, one 
solves the Kohn-Sham equations using the approxima¬ 
tion of Eq. ([5| in order to obtain a new density. The pro- 
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FIG. 1: (color online) Densities for N = 4 dipoles in a Q1D 
geometry ( l x = 101x) as a function of the dipole moment 
d , which is assumed to be oriented perpendicularly to the 
direction of motion. The different results correspond to the 
KS-SCP DFT (red for bosons (B), black for fermions (F)), 
Cl (blue for bosons, green for fermions), and CP (magenta) 
approaches. 

cess is then repeated self-consistently until convergence 
is achieved. 

To illustrate the kind of accuracy that can be obtained 
with the KS-SCP formalism, we chose a case for which 
full Cl calculations can be also performed: in Fig.[l]we re¬ 
port the KS-SCP densities for a system with four bosonic 
and fermionic dipoles in a Q1D trap. To model this sys¬ 
tem we have chosen the oscillator lengths in the longitu¬ 
dinal and perpendicular directions l x = yjh/ ( muj x ) and 
l± = i Jh/ (mwi) such that l x = 10Z_l- Also, we have 
assumed the dipoles to be oriented perpendicularly to 
their direction of motion, that is 0 r d = 0 = 90°. For 
comparison, we also show the results obtained from the 
mean-field GP approach for the bosonic case. When d 
is very small (0.01 and 0.1) the system is in the weakly- 
interacting regime: we see good qualitative agreement 
of the KS-SCP results with the Cl for fermions, and 
also with the GP calculations for bosons, reflecting the 
fact that the correlation effects are negligible. As d is 
further increased, however, the picture changes qualita¬ 
tively. For d = 1 the Cl bosonic results show a den¬ 
sity structure characteristic of the so-called fermionized 
regime [2j3], with two tall central peaks and two shorter 
lateral ones, indicating that the bosons feel the infinitely 
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strong short-ranged part of the interaction and behave 
as non-interacting spinless fermions. The zeroth-order 
KS-SCP approach can only describe this phenomenon in 
a qualitative way, barely displaying two incipient lateral 
peaks but being not able to reproduce the central struc¬ 
ture for the chosen parameters. This is due to the large 
underestimation of the kinetic correlation energy: in the 
SCP state of Eq. Q, the particles are located in their 
strictly-correlated positions, minimizing the interaction 
energy for the given density, but without increasing sub¬ 
stantially the kinetic energy, which is obtained from the 
KS construction. The resulting SCP potential of Eq. Q 
is then too small, since the particle are too far from each 
other. Still, one can see that the KS-SCP results are 
much better than those obtained from the GP approach: 
indeed, the latter yields a Thomas-Fermi-like density pro¬ 
file lacking of any structure. Finally, for strong enough 
values of the dipole moment (d = 8, 20 and 50), the 
system enters the localized regime, where the Cl densi¬ 
ties show a characteristic profile with four clearly marked 
peaks corresponding to the localization of the density [25] 
due to the strong long-range repulsion between the par¬ 
ticles. One can see how the KS-SCP densities show this 
structure as well, and that they become closer to the ex¬ 
act ones as the strength of the interaction increases. The 
capacity of the KS SCP approach for going beyond the 
mean-field description is apparent from the case d = 8, 
where the GP density illustrates the total neglect of cor¬ 
relation effects (for larger d we could not even get con¬ 
verged GP results within our grid). It is also remarkable 
that the same approximate method is able to span (even 
if only qualitatively) a wide range of different correlation 
regimes. Moreover, as any DFT approach, the formalism 
is amenable to corrections: exchange effects, for exam¬ 
ple, could be included in the SCP functional in order to 
capture the fermionized regime. Notice that while full 
Cl calculations almost reached the maximum number of 
particles that can be treated, the computational cost of 
the KS-SCP method in one dimension is similar to that 
of the usual LDA. 

In Fig. [2] we show the KS-SCP densities for Q2D sys¬ 
tems with N = 15 particles, considering the dipole mo¬ 
ments to be perpendicularly aligned with the plane of 
motion, and setting 10l z = l± = l x = l y = 1. Sim¬ 
ilarly as in the Q1D case of Fig. |T| one can see how 
the KS SCP method is able to span the range from 
the weakly-interacting regime (corresponding to Bose- 
Einstein condensation for bosons and to the shell struc¬ 
ture for fermions), to the strongly-correlated one. In this 
latter case, one can observe the formation of two char¬ 
acteristic concentric rings, which by integration of the 
density reveal the presence of 6 and 9 particles on aver¬ 
age. Again, these calculations are out of reach of full Cl, 
and also difficult with QMC. 

Conclusions - In this Letter we have suggested an 
exchange-correlation functional for the application of 
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FIG. 2: (color online) KS-SCP densities for bosonic and 
fermionic systems in a Q2D geometry for N = 15 in the 
weakly- and strongly-correlated regimes. 

Kohn-Sham DFT to ultracold bosonic and fermionic 
dipolar and ionic gases. This functional is based on 
the exact strong-coupling limit for a given density n(r), 
and can be used to treat long-ranged interactions in the 
strongly-correlated regime. The functional can be im¬ 
proved with rigorous corrections (33J [56] and one should 
expect that inclusion of higher-order terms would im¬ 
prove the performance for intermediate regimes. 

The results show the potential of strictly-correlated- 
based density functionals for the study of ultracold dipo¬ 
lar and ionic gases, where only little progress has been 
done in the application of density-functional methods m 
compared to ultracold atomic systems with short-ranged 
interactions. The unique properties and high tunability 
of the dipole-dipole interaction have been shown to be po¬ 
tentially useful for the study of many interesting phenom¬ 
ena, as well as for several practical applications. Their in¬ 
vestigation with density-functional methods might open 
new possibilities, especially when generalized to the time 
domain (which is the object of ongoing work). The pre¬ 
sented method can be applied to systems with other ex¬ 
ternal potentials, such as the interesting case of ID opti¬ 
cal lattices, as well as generalized to different interactions 
and higher dimensionalities. In this sense, one interesting 
possibility could be the recently discussed quadrupolar 
interactions, which may lead to intriguing new quantum 
phases (see [5%HSH1 ). These topics will be the subject of 
future works. 

Acknowledgments - We thank D. Karlsson, 
G. Kavoulakis, I. Tokatly, C. Verdozzi and G. Vig- 
nale for fruitful discussions. We acknowledge support 
from the Netherlands Organization for Scientific Re¬ 
search (NWO) through a Vidi grant, from a Marie Curie 
Intra European Fellowship within the 7th European 




5 


Community Framework Programme, from the Swedish 
Research Council, and from the Nanometer Structure 
Consortium at Lund University. 


[1] A. J. Cohen, P. Mori-Sanchez, and W. Yang, Chem. Rev. 
112, 289 (2012). 

[2] J. Neugebauer and T. Hickel, Wiley Interdisciplinary Re¬ 
views: Computational Molecular Science 3, 438 (2013). 

[3] S. M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002). 

[4] P. Hohenberg and W. Kohn, Phys. Rev. 136, B 864 
(1964). 

[5] W. Kohn and L. J. Sham, Phys. Rev. A 140, 1133 (1965). 
6 I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 

[7] A. Minguzzi, S. Sued, F. Toschi, M. Tosi, and P. Vignolo, 
Phys. Rep. 395, 223 (2004). 

[8] L. P. P. S. Giorgini and S. Stringari, Rev. Mod. Phys. 
80, 1215 (2008). 

[9] M. Rontani, C. Cavazzoni, D. Bellucci, and G. Goldoni, 
J. Chem. Phys. 124, 124102 (2006). 

[10] E. Waltersson, C. J. Wesslen, and E. Lindroth, Phys. 
Rev. B 87, 035112 (2013). 

[11] A. Ghosal, A. D. Guclu, C. J. Umrigar, D. Ullmo, and 
H. U. Baranger, Nature Phys. 2, 336 (2006). 

[12] A. Ghosal, A. D. Guclu, C. J. Umrigar, D. Ullmo, and 
H. U. Baranger, Phys. Rev. B 76, 085341 (2007). 

[13] G. S. Nunes, J. Phys. B: At. Mol. Opt. Phys. 32, 4293 
(1999). 

[14] Y. E. Kim and A. L. Zubarev, Phys. Rev. A 67, 015602 
(2003). 

[15] J. Brand, J. Phys. B: At. Mol. Opt. Phys. 37, 287 (2004). 

[16] Y. Hao and S. Chen, Phys. Rev. A 80, 043608 (2009). 

[17] P. N. Ma, S. Pilati, M. Troyer, and X. Dai, Nature Phys. 
8, 601 (2012). 

[18] M. Casula, S. Sorella, and G. Senatore, Phys. Rev. B 74, 
245427 (2006). 

[19] N. A. Lima, M. F. Silva, L. N. Oliveira, and K. Capelle, 
Phys. Rev. Lett. 90, 146402 (2003). 

[20] H. Wang, Y. Hao, and Y. Zhang, Phys. Rev. A 85, 053630 

( 2012 ). 

[21] H. Wang and Y. Zhang, Phys. Rev. A 88, 023626 (2013). 

[22] M. A. Baranov, Physics Reports 464, 71 (2008). 

[23] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and 
T. Pfau, Rep. Prog. Phys. 72, 126401 (2009). 

[24] M. A. Baranov, M. Dalmonte, G. Pupillo, and P. Zoller, 
Chem. Rev. 112 , 5012 (2012). 

[25] K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe’er, 
B. Neyenhuis, , .1. J. Zirbel, S. Kotochigova, P. S. Juli¬ 
enne, D. S. Jinl, et ah, Science 322, 231 (2008). 

[26] S. Ospelkaus, K.-K. Ni, D. Wang, M. H. G. de Miranda, 
B. Neyenhuis, G. Quemener, P. S. Julienne, J. L. Bohn, 
D. S. Jin, and .J. Ye, Science 327, 853 (2010). 

[27] G. E. Astrakharc.hik, G. Morigi, G. DeChiara, and 
J. Boronat, Pliys. Rev. A 78, 063622 (2008). 

[28] N. T. Zinner, B. Wunsch, I. Mekhov, S.-J. Huang, 
D. Wang, and E. Dernier, Phys. Rev. A 84, 063606 
( 2011 ). 

[29] F. Deuretzbacher, J. C. Cremon, and S. M. Reimann, 


Phys. Rev. A. 81, 063616 (2010). 

[30] M. Dalmonte, G. Pupillo, and P. Zoller, Phys. Rev. Lett. 
105, 140401 (2010). 

[31] R. Citro, E. Orignac, S. D. Palo, and M. Chiofalo, Phys. 
Rev. A 75, 051602(R) (2007). 

[32] B. Wunsch, N. Zinner, I. Mekov, S.-J. Huang, D.-W. 
Wang, and E. Demler, Phys. Rev. Lett. 107, 073201 
( 2011 ). 

[33[ P. Gori-Giorgi, M. Seidl, and G. Vignale, Phys. Rev. Lett. 
103, 166402 (2009). 

[34] P. Gori-Giorgi and M. Seidl, Phys. Chem. Chem. Phys. 
12, 14405 (2010). 

[35] F. Malet and P. Gori-Giorgi, Phys. Rev. Lett. 109, 
246402 (2012). 

[36] F. Malet, A. Mirtschink, J. C. Cremon, S. M. Reimann, 
and P. Gori-Giorgi, Phys. Rev. B 87, 115146 (2013). 

]37] C. B. Mendl, F. Malet, and P. Gori-Giorgi, Phys. Rev. B 
89, 125106 (2014). 

[38] E. H. Lieb, Int. J. Quantum. Chem. 24, 24 (1983). 

[39] M. Seidl, Phys. Rev. A 60, 4387 (1999). 

[40] M. Seidl, P. Gori-Giorgi, and A. Savin, Phys. Rev. A 75, 
042511 (2007). 

[41] G. Buttazzo, L. Dc Pascale, and P. Gori-Giorgi, Phys. 
Rev. A 85, 062502 (2012). 

[42] C. Cotar, G. Friesecke, and C. Kliippelberg, Comm. Pure 
Appl. Math. 66, 548 (2013). 

[43] C. Cotar, G. Friesecke, and B. Pass, in Calculus of Varia¬ 
tions and Partial Differential Equations (Springer, Berlin 
Heidelberg, 2014), pp. 1-26. 

[44] G. Friesecke, C. B. Mendl, B. Pass, C. Cotar, and 
C. Kliippelberg, J. Chem. Phys. 139, 164109 (2013). 

[45] M. Colombo, L. De Pascale, and S. Di Marino, Can. J. 
Math. , accepted (2014). 

[46] B. Pass, Nonlinearity 26, 2731 (2013). 

[47] B. Pass, J. Func. Analysis 264, 947 (2013). 

[48] N. Ghoussoub and A. Moameni, Geom. Funct. Anal. 24, 
1129 (2014). 

[49] Y.-H. Kim and B. Pass, SIAM J. Math. Anal. 46, 1538 
(2014). 

[50] M. Colombo and S. Di Marino, in Annali di Matematica 
Pura ad Applicata (Springer, Berlin Heidelberg, 2013), 
pp. 1 14. 

[51] A. Mirtschink, M. Seidl, and P. Gori-Giorgi, J. Chem. 
Theory Comput. 8, 3097 (2012). 

[52] E. Rasanen, M. Seidl, and P. Gori-Giorgi, Phys. Rev. B 
83, 195111 (2011). 

[53] J. C. Cremon, G. M. Bruun, and S. M. Reimann, Phys. 
Rev. Lett. 105, 255301 (2010). 

[54] C. B. Mendl and L. Lin, Phys. Rev. B 87, 125106 (2013). 

[55] H. Chen, G. Friesecke, and C. B. Mendl, J. Chem. Theory 
Comput. 10 , 4360 (2014). 

[56] P. Gori-Giorgi, G. Vignale, and M. Seidl, J. Chem. The¬ 
ory Comput. 5, 743 (2009). 

[57] B. Fang and B.-G. Englert, Phys. Rev. A 83, 052517 

( 2011 ). 

[58] M. Lahrz, M. Lemeshko, and L. Mathey, 

arXiv:1412.8472S (2014). 

[59] S. G. Bhongale, L. Mathey, E. Zhao, S. Yelin, and 
M. Lemeshko, Phys. Rev. Lett. 110 , 155301 (2013). 

[60] W.-M. Huang, M. Lahrz, and L. Mathey, 

arXiv:1311.1947v2 (2014). 



